Phase diagram of the one-dimensional, two-channel Kondo lattice model 
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Employing the density matrix renormalization group method and strong-coupling perturbation 
theory, we study the phase diagram of the SU(2) xSU(2) Kondo lattice model in one dimension. We 
show that, at quarter filling, the system can exist in two phases depending on the coupling strength. 
The weak-coupling phase is dominated by RKKY exchange correlations while the strong-coupling 
phase is characterized by strong antiferromagnetic correlations of the channel degree of freedom. 
These two phases are separated by a quantum critical point. For conduction-band fillings of less 
than one quarter, we find a paramagnetic metallic phase at weak coupling and a ferromagnetic phase 
at moderate to strong coupling. 
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Although Landau's theory of Fermi hquids is one of 
the cornerstones of modern condensed matter physics," 
many materials show metallic properties that do not fit 
into this framework. For instance, it is well known that 
one-dimensional systems usually behave as Luttinger liq- 
uids rather than Fermi liquids.^ In higher dimensions, the 
proximity of a quantum critical point (QCP) is consid- 
ered to be responsible for the non-Fermi liquid properties 
of systems like CeCue-aAua;!^ The quadrupolar Kondo 
effect has been proposed as an alternative source of non- 
Fermi liquid behavior. This effect is described by the two- 
channel Kondo model^. This model has a non-Fermi 
liquid ground-state because of frustration in the screen- 
ing of a localized impurity by two degenerate conduction 
electron channels if the degeneracy of the channels M 
is greater than twice the impurity spin S.^ The single- 
impurity two-channel Kondo model can provide an ade- 
quate description of dilute systems like Thi_a;Ua;Ru2Si2^ 
or Yi_a;Ua;Pd3(S but it does not incorporate the lattice 
effects that become relevant in fully concentrated com- 
pounds like UBcis^. These materials are described by 
the two-channel Kondo lattice model (KLM). 

In this paper, we address the question of what happens 
when two of these fundamental sources of non- Fermi liq- 
uid behavior coincide. For this purpose we study the 
SU(2)xSU(2) Kondo lattice model in one spatial dimen- 
sion. The Hamiltonian reads 



H = -t 
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where t > is the conduction electron hopping ampli- 
tude, taken to be the same in both bands, and cj„j^ {Cima) 
creates (annihilates) an electron on lattice site i, I < i < 
L {L being the number of lattice sites), with channel fla- 
vor m = + or — and spin projection a =t or |. The 
Heisenberg spin operator for the localized /-electrons 



i) is defined by = 5 //a'^a/3/»/3' where 



the /-operators satisfy the constraint //i/^i + //|/i| = 1 
and cr is a vector of Pauli spin matrices. The conduction 
band filling is defined as rif. = n^+^n^.^ = [Nc++Nc-) / L 
with Nc± the number of conduction electrons in channel 
m — zL, respectively; ric — 1 corresponds to the quar- 
ter filled system. The KLM may be derived from the 
more fundamental periodic Anderson model in the limit 
of strong Coulomb repulsion where the Kondo coupling 
is usually antiferromagnetic (AF), J > 0,^ or alterna- 
tively in the "extended Kondo limit" In the following, 
we measure all energies in units of t. 

Not much is known about the ground-state phase di- 
agram of the two-channel KLM. Tsvelik and Venturaii 
investigated this model using a mean-field analysis and 
found that at half filling the system exists in two phases. 
One is dominated by RKKY exchange interaction ef- 
fects, and the other by Kondo screening. A QCP sep- 
arates these two regimes. A generalized one-dimensional 
two-channel KLM with an additional Heisenberg interac- 
tion, Jh, between the /-spins was studied by Andrei and 
Orignac."^ In the limit of strong Jh, they find that the 
system is in a superconducting phase with odd-frequency 
singlet pairing of the electrons. In infinite spatial dimen- 
sions, the ground state is characterized by superconduct- 
ing or magnetic phases, which may coexist or competed 

The focus on the iwo-channel system below or at quar- 
ter filling is also motivated by the relation to the single- 
channel system at half filling or less. In particular, the 
quarter-filled case for the two-channel model is analogous 
to the half-filled case for the single-channel model in that 
there is one conduction electron per impurity spin, which 
leads to complete screening at strong coupling. In the fol- 
lowing, we will show that the situation in the two-channel 
model below quarter filling is qualitatively similar to the 
single-channel case below half filling, while exactly at 
quarter filling it is quite different. The phase diagram 
of the single-channel model is well understood, at least 
qualitativelyi^*i^ In the low carrier limit, this system 
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displays ferromagnetic order with complete polarization, 
5tot = (i - Nc)/2, where < In the strong- 

coupling limit, the ground state is ferromagnetic for all 
rici^ Exactly at half filling, the single-channel model is 
known to be a Kondo insulator. ^'^•^^ This is a quantum 
disordered phase in which the conduction electrons are 
bound into local singlet states with the impurity spins, 
and both the spin and charge correlation functions decay 
exponentially in space and time. 

In order to calculate the ground-state properties of 
the one-dimensional two-channel KLM, we use the finite- 
system algorithms^ of the DMRG to calculate gaps, 
equal-time correlation functions, and the total spin of 
the ground state. We keep up to 1000 states per block 
on lattices of up to L = 50 sites and obtain a maximum 
discarded weight of 10~^. Fig. shows the total spin 
'S'tot per site extrapolated to the thermodynamic limit 
for various couplings and conduction electron densities. 
Here 5'tot is calculated directly by taking the expecta- 
tion value of Sjoj in the ground state and also by ex- 
amining the degeneracy of excited states in various Sz 
sectors. A finite-size extrapolation is then used to de- 
termine whether Stot = {L — Nc)/2 (complete ferromag- 
netism), S'tot < {L — Nc)/2 but finite (incomplete ferro- 
magnetism) , or S'tot — (paramagnetism) in the thermo- 
dynamic limit. At quarter filling (ric = l)j the nature of 
the ordering is also indicated for 5tot — phases. 
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FIG. 1: Ground-state phase diagram of the two-channel KLM 
as a function of conduction band filling Uc- The AF channel- 
and spin-ordered phases at quarter filling are associated with 
a singlet ground state. The crosses indicate polarizations that 
extrapolate in the thermodynamic limit to values between 
25% and 90% of complete polarization. 

Fig-fflshows a large region of complete ferromagnetism 
for He < 1 above a certain critical value Jc which de- 
creases with decreasing Uc and tends to zero as nc ^ 0. 
In the completely polarized phase, each conduction elec- 
tron forms an itinerant singlet with the /-spins. Delocal- 
ization of these singlets leads to ferromagnetic ordering 
of the remaining unscreened /-spins. The complete po- 
larization for all J at low conduction electron density is 
in agreement with an exact argument for a single con- 



duction electron. For J < Jc, there is a region of incom- 
plete ferromagnetism, similar to one that has been found 
in the periodic Anderson modeliSl While we cannot rule 
out that this region is due to a continuous transition to 
the complete ferromagnetic phase, the local spin profiles 
in the incomplete ferromagnetic regime show small fer- 
romagnetic domains (corresponding to polarizations be- 
tween 25% and 90% of the complete value), suggesting 
that phase separation may occur here. For ric > 1, we 
find a singlet ground state for all couplings J. 

Exactly at quarter filling, we find two phases as a func- 
tion of J. At weak coupling (J < 2.0), the electrons 
of different flavors generate independent RKKY interac- 
tions between the localized moments. We observe strong 
correlations of the /-spins at a wave vector q — tt/2, as in 
the single-channel model at quarter filling. At stronger 
coupling ( J 2.0), the system is in a channel AF phase, 
where the correlations of the channel degree of freedom 
decay as 1/r. The two phases are separated by a QCP. 
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FIG. 2: Magnetic correlation function S/(r) = {5'z(0)S'z(r)) 
of the /-spins for the two-channel KLM at tIc = 1 in (a) and 
for the single- and two-channel models, both at quarter-filling, 
in (b). 

Fig.[2Ia) shows that in the two-channel model the spin- 
spin correlation function decays slowly for weak coupling 
(J = 0.5 or 1.0) and is short-ranged for strong couphng 
(J = 4.0), i.e., numerically zero for more than two lat- 
tice spacings. The transition between the two behaviors 
occurs at J « 2.0 (not shown in Fig. where the cor- 
relations extend over roughly two lattice sites. Fig. |2b) 
shows a comparison between the spin-spin correlations 
functions Sf{r) of the single-channel and two-channel 
KLM on a logarithmic plot. The ground state for M = 1 
in the weak-coupling limit shows RKKY liquid behavior 
with spatial oscillation characterized by a wave vector 
q = tt/2. In the M = 2 case the correlations at J = 0.5 
decay so slowly that the asymptotic behavior cannot be 
determined for the system sizes considered, as can be 
seen in Fig. |2Ib). At larger J- values [e.g., J = 1.0 in 
Fig. OJb)] , the correlation functions appear to decay ex- 
ponentially, consistent with the opening of a gap in the 
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spin excitation spectrum due to the transition into the 
channel-ordered phase. Fig. |21 shows the magnetic struc- 
ture factor Sf{q) for the two-channel system. The am- 
plitude of the peak a,t q — tt /2 (the wavevector expected 
for a RKKY liquid) becomes smaller with increasing J, 
indicating that the magnetic correlations between the /- 
spins become weaker and finally vanish at J ~ 2.0m^ 

It is known that the RKKY correlations of the 
single-channel KLM at quarter filling become unstable 
with increasing coupling toward a ferromagnetic ground 
statejiS*^ However, the two-channel model cannot be 
ferromagnetic at quarter filling because each /-spin is 
screened by an electron, so that some other kind of sym- 
metry breaking must lift the degeneracy of the ground 
state. Fig. 21 shows the staggered magnetization D{r) of 
the channel degree of freedom, with 
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4,naC„na- In Fig- Ufa) one sees that D{r) 



where riima 

decays with 1/r as function of the distance r — i ~ j for 
J = 4.0. Fig.^Jb) shows that the AF channel correlations 
become stronger with increasing coupling and the quasi- 
long-range behavior develops for J > 2.0. 
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FIG. 3: Magnetic structure factor of the two-channel KLM 
model as a function of wavevector q for various J at ric = 1.0. 

In order to understand our numerical results for strong 
coupling, we use the methods of Ref. to derive an 
effective Hamiltonian valid for i/J <C 1. At ric = 1, 
the channel flavor is the only degree of freedom and the 
low-energy spectrum can be described by a pseudospin- 
1/2 model. Accordingly, the effective Hamiltonian is a 
Heisenberg model for the channel spin. 
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where r = \ Y.r,im' ^iTn^rn-m'b^m' defined in terms of 
hard-core bosons 



(2) 



which represent local Kondo singlets. In one spatial di- 
mension, model Q has been solved exactly using the 
Bethe ansatzi^^ The ground state is characterized by crit- 
ical AF correlations that decay as 1/r, in agreement with 
our numerical findings at quarter filling [see Fig. 2Ja)]. 
This critical behavior is replaced by long-range AF or- 
dering for D > 1. Due to the SU(2) pseudo-spin in- 
variance of H, the AF ordering is present for all three 
r-components. In particular, the order along the z axis 
corresponds to staggered orbital ordering while the x, y- 
ordering gives rise to a Bose-Einstein condensation of ex- 
citons (particle-hole singlet bound states between the two 
channels). 

For Tic < 1, the low-energy subspace in the strong- 
coupling limit contains states with zero and one conduc- 
tion electron per site. While the local state with one 
conduction electron can be described by the hard-core 
bosons of Eq. the empty state acquires spin-1/2 
character from the / spins. It can therefore be rep- 
resented as a constrained fermion (p]^ = (1 — n'^)f^„, 

with ~ J2mbimbim- At strong couphng, the effective 
Hamiltonian has the form 



Ho 



i(7,m 



In analogy to the infinite-?/ Hubbard model, the ground- 
state wavefunction can be written as a direct product of 
a charge, spin, and orbital component 



\'4'n{T; a}) ^ \nc) In - - - tl^n) ®\(ti - ■■ (tn) 
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where {t; cr} = (ti, r^.^r; Ci, cttv) and |0) denotes 
the vacuum of ip and b particles. The complete spin de- 
generacy of -ffo in the strong coupling limit is lifted in 
0{t^/J). The effective Hamiltonian in this order con- 
tains only one term that lifts this degeneracy: 



t / 



where nf = X^o- '^L'^io-- Here Hi exchanges two spins 
by hopping of a fermion over another to an empty next- 
nearest-neighbor site. In the same way as for the single- 
channel KLM, it can be shown that the off-diagonal ele- 
ments of Hi are all non-positive and the Hamiltonian ma- 
trix in real-space is completely connected. The Perron- 
Frobenius theorem (see, for instance, Ref. l2^ then states 
that the ground state is unique and that the coefficients 
of the wave-function can be chosen to be strictly positive. 
The only spin state that has strictly positive coefficients 
in each of the subspaces of Hi is the one with maximum 
total spin. Therefore, the ground state of the two-channel 
KLM below quarter-filling is ferromagnetic in the strong- 
coupling limit, in accordance with our numerical findings. 
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FIG. 4: Staggered magnetization D{r) for the channel degree 
of freedom for various J at quarter filling {ric = 1.0). The 
longer distance behavior for J = 4.0 (where there is a slower 
falloff, approximately oc 1/r) is shown in (a) and a shorter 
range in r is shown for all four J- values in (b). 



In this paper, we have mapped out the zero-tempera- 
ture phase diagram of the two-channel KLM as a func- 
tion of conduction band filling ric and Kondo coupling 
strength J. Our main results are that the phase dia- 
grams of the two-channel and single-channel KLMs are 



qualitatively similar for low band fillings (ric < 1 in both 
models) but quite different at = 1. For < 1, the 
addition of a second degenerate band of conduction elec- 
trons does not alter the physical picture of the single- 
channel model if ric is low or J is large. The similarity 
between the two phase diagrams for ric < 1 further sug- 
gests that the metallic ground state of the two-channel 
KLM at weak coupling is determined by two independent 
Kondo effects in both channels. In contrast, at quarter- 
filling (ric = 1) the two-channel KLM exhibits a quantum 
phase transition between this metallic phase and an in- 
sulator characterized by strong AF correlations of the 
channel degrees of freedom. Our perturbative analysis 
shows that this insulating phase is present for any spa- 
tial dimension D and has long-range ordering for D > 1. 
Therefore, we also expect a QCP separating the metal- 
lic and the insulating phase for D > 1. The AF chan- 
nel (or excitonic) fluctuations diverge at the QCP and 
can produce a deviation from the normal Fermi-liquid 
behavior of the metallic state. In this way we see that 
the deviations from the Fermi-liquid behavior which are 
obtained in the dilute (impurity model) and the concen- 
trated (lattice model) limits have a common origin in the 
fluctuations of the channel degree of freedom. 
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